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Disoriented Chiral Condensates (DCC) have been predicted to form in high energy 
heavy ion collisions where the approximate chiral symmetry of QCD has been restored. 
This leads to large imbalances in the production of charged to neutral pions. Sophisti- 
cated analysis methods are being developed to disentangle DCC events out of the large 
background of events with conventionally produced particles. We present a short review 
of current analysis methods and future prospects. 



1. Introduction 

The QCD phase transition from normal hadronic matter to the Quark-Gluon-Plasma 
(QGP), in case of high energy heavy ion collisions, manifests itself in two forms: (1) 
Deconfmement transition and (2) Chiral symmetry restoration. One of the interesting 
consequences of the chiral transition is the formation of a chiral condensate in an extended 
domain, such that the direction of the condensate is misaligned from the true vacuum 
direction. This is called the formation of Disoriented Chiral Condensates (DCC). It has 
been proposed that the decay of the DCC domains would lead to large imbalances in 
the production of charged to neutral pions. The task for experimentalists is to carefully 
measure the number of neutral and charged pions as well as study their spectra. The 
challenge is to design sophisticated analysis tools on an event-by-event basis to identify 
DCC amidst the large background due to conventionally produced particles. 

The formation of DCC domains has been proposed by Anselm[ ffl], by Blaizot and 
Krzywicki[ |2] and by Bjorken, Kowalski and Taylor [ |3| in the context of high energy 
hadronic collisions in order to explain the puzzling Centauro and anti-Centauro type of 
events observed by cosmic ray experimentsf Bjorken et al.[ ^ [6[] proposed the so called 
"Baked Alaska" model which suggests that in case of hadronic collisions a hot fireball is 
produced with cold interiors having anomalous chiral order parameter. Rajagopal and 
Wilczek[ U were the first ones to discuss the DCC phenomena in the context of heavy 
ion reactions. They have suggested that the nonequilibrium dynamics during the chiral 
symmetry breaking phase transitions in case of heavy ion collisions may produce DCC 
domains. There is tremendous progress [ @, 0, £| H [H], Q; * n t erms °f theoretical 



understanding of DCC since then starting with different formalisms. 

In the framework of the linear sigma model the Lagrangian can be expressed in terms 
of the order parameter $ = (a, 7?) , which is a combination of the scalar field a and the 
pion field ff. At low temperatures the chiral symmetry is spontaneously broken. The 
potential, V(<&), has a minimum in the sigma direction and all the direction of tt are 
equally populated. Thus the distribution of neutral pion fraction defined by 

/ = — (1) 

is a gaussian with a mean of 1/3. This leads to the isospin symmetry of pions. As the 
system goes through the chiral transition and then rapidly expands and cools, it may roll 
down from the unstable local maximum of V(<&) to the nearly stable values of one of the 
pion directions. The field will have to eventually settle in to the true ground state, but 
oscillations will continue for sometime which leads to amplification of soft pion modes. 
This effectively creates a cluster of low pt pions in a correlated region (of a so called DCC 
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domain) with the probability of the neutral pion fraction given by: 

P(f) - (2) 

As is evident from this discussion, pions from a DCC domain will be emitted at low 
Pt and will have a distinct distribution pattern compared to the normal pion production 
mechanism without DCC. 

Our ability to detect DCC domains depends on the number of domains, size of domains 
and number of DCC pions emitted from the domains. If the number of domains is more 
than one, then the distribution given by (2) gets modified, and for a large number of 
uncorrelated small domains, the resulting distribution becomes gaussian. The size of a 
DCC domain and the energy content determine the number of DCC pions. Clearly the 
fewer the number of domains and larger the number of emitted pions, the easier it is 
to detect in the laboratory. In addition, the probability of DCC formation in a nuclear 
collision at a given energy is important for the observation of DCC. 

The cleanest evidence for the DCC formation can be obtained by precision measurement 
of N w o and N n ±, or equivalently, iV 7 and distributions. This assumes all charged 
particles are pions and all photons come from tt° decays. In this manuscript we will 
deal mostly with analysis of neutral and charged particle distributions. Other signals 
which have been proposed for detection of DCC will be discussed towards the end of the 
manuscript. 



2. Accelerator Experiments 

The observation of centauro type of events have attracted lot of attention since the 
early 80's. Several accelerator experiments have been performed to search for such un- 
usual events even before the concept of "DCC" came into being. The UAl and UA5 
experiments at the CERN-SPS have carried out search for centauros in p — p collisions 
at CM energies of 540 and 900 GeV, respectively, whereas DO and CDF experiments at 
Fermilab have performed this search for CM energy of 1.8 TeV. These experiments have 
used the technique of the asymmetry of hadronic to electromagnetic energies. So far there 
is no evidence of any centauro type of events from these experiments. 

Below we describe two current experiments for DCC search: the Minimax experiment 
for p — p collisions and the WA98 experiment for heavy ion collisions. Results from WA98 
experiment will be discussed in detail. Other experiments at CERN-SPS such as NA45 
and NA49 have plans for DCC search in near future. 

2.1. Minimax Experiment at Fermilab 

The Minimax Experiment [ O] at Fermilab Tevatron is designed by Bjorken et al. for 
the DCC search in p—p collisions of 1.8TeV in the far forward region. The experiment is 
located at the CO interaction region of the Tevatron and has been designed to measure 
the ratio of charged to neutral pions produced at pseudo-rapidities near 4.1. The lego 
acceptance of the detector is a circle of radius 0.65 units in rj. 

The data analysis for DCC search in Minimax is complicated because of small accep- 
tance of the detector and various efficiency factors which come into play. The method 
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of normalized factorial moments has been utilized successfully to construct a set of "ro- 
bust observables" which eliminate most of these effects [ [13| . This method is described in 
section 5.3. The results from these data is consistent with no DCC production mechanism. 

2.2. WA98 Experiment at CERN-SPS 

The WA98 experiment [ |14| at CERN-SPS emphasizes on high precision, simultaneous 
measurement of both hadrons and photons. The experimental setup consists of large 
acceptance hadron and photon spectrometers, detectors for charged particle and photon 
multiplicity measurements, and calorimeters for transverse and forward energy measure- 
ments. At present our search is limited to detailed event-by-event analysis of photon 
distributions from the Photon Multiplicity Detector and charged particles from Silicon 
Pad Multiplicity Detectors. Details of these detectors may be found elsewhere, here we 
give the essential points necessary for our present discussion. 

The Photon Multiplicity Detector (PMD) 

The PMD consists of a 3X thick lead converter sheet in front of an array of 54,000 
pads, and is located at a distance of 21.5 meters from the target. Signals from several 
neighboring pads are combined to form clusters, characterized by the total ADC content 
and the hit positions. A threshold of 3 MIPs on the preshower ADC content gives an 
average photon counting efficiency in the range of 70%-75% depending on the centrality, 
with about 32% contamination due to showering hadrons. Thus the measured particle 
clusters, called "7-like clusters" contain mostly photons together with a sizable contribu- 
tion of charged pions. 

The Silicon Pad Multiplicity Detector (SPMD) 

The SPMD is located at 32.8 cm from the target. It is based on a double metal, AC 
coupled design consisting of approximately 4000 pads, arranged radially with 180 0— bins 
and 22 equal— rj bins between 7] = 2.35 and 3.75. The granularity of the detector is high 
enough that we have only 15 — 20% occupancy in the highest multiplicity events. The 
efficiency of detecting a charged particle in the active area has been determined in a test 
beam to be 99%. Conversely the detector is transparent to high energy photons, since 
only about 0.2% are expected to interact in the silicon. 

3. DCC Analysis Methods 

Most of the DCC analyses so far are based on studying the fluctuations in neutral 
and charged particle measurements. Various analysis methods are being employed to 
characterize unusual events which show up beyond the statistical fluctuations. All the 
analysis in WA98 are being carried out in event-by-event basis. Current techniques in 
DCC analysis are the following: 

• Global event characteristics: This analysis is performed by using the total number 
of photons and charged particles over the entire phase space covered by the photon 
and charged particle detectors. 

• Methods for DCC domains localized in phase space: The available phase space is 
divided into several rj — <fi bins. Here we mention two of the analysis methods which 
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are sensitive enough to single out localized domains: 

- Wavelet Analysis: This multiresolution analysis has the capability to scan the 
entire phase space. No averaging over events or r\ — <p space is done. 

- Moments analysis: Various moments and combinations of those are calculated 
from the distribution of photons and charged particles in each bin. 

In all these analysis, data are not corrected for detector effects and other efficiencies. 
Instead we compare the data with MC simulations which incorporates all known detector 
and physics effects. The output of VENUS 4.12 event generator is passed through a full 
simulation of our experimental setup using the GEANT 3.21 package from CERN. The 
same method of analysis is followed for both data and MC generated events. The idea is 
to look for 

• shapes of different distributions: non-gaussian shapes in case where gaussian is 
expected and differences in various distributions in case of data and MC simulation 

• exotic events: in case the shape is deviating from what is expected, one can study 
the events which show up beyond some predetermined limit. 

4. Global Event Characteristics 

The method of global event characterization in terms of the photon and charged particle 
multiplicity distributions over the full available phase space is suitable for the search of 
single large size DCC domain. The idea here is to look for events which fall far beyond 
the correlation line of these two distributions. For the current study we have assumed 
that the DCC could be formed only in central collisions. We will concentrate on the 
10% most central events, defined by a measured transverse energy of at least 300 GeV in 
3.5 < rj < 5.5. 

In Figures la and lb we present the minimum-bias multiplicity distribution for 7-like 
clusters and charged particles. The central data sample is shown by filled circles in 
these two figures. After all cuts are applied, there are 212646 events in this sample. A 
comparison with the MC simulation events chosen by identical cuts for central events is 
shown by the histogram. 

The correlation between the charged and neutral multiplicities is presented on the right 
side of Figure 1 with the minimum bias distribution outlined, the central MC simulated 
events hatched, and the central data events shown as scattered points, each point cor- 
responding to a single event. A strong correlation is seen between charged and neutral 
multiplicities, which suggests a more appropriate coordinate system with one axis being 
the measured correlation axis and the other perpendicular to it. This is represented by 
the Z axis and the Dz axis as shown in figure lb. We have chosen to work with the scaled 
variable Sz=Dz/cd z i n order to compare relative fluctuations at different multiplicities. 
Sz distribution for the data is shown as filled circles in figure 2(a). 

4.1. DCC simulation 

To estimate the effect of DCC production we have modified the output of VENUS 
events to include characteristic fluctuations in the relative production of charged and 
neutral pions. We assume that only a single domain of DCC is formed in each central 
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Figure 1. Multiplicity distributions of (a) 7-like clusters and (b) charged particles for 
minimum-bias (open circles) and central events (filled circles). The histograms show MC 
simulation results for central events. The right hand panel shows the correlation of iV 7 _iik c 
and N c h for central events. The hatched region is for MC simulation results. 



collision. A certain fraction, £, of the VENUS pions is associated with this domain and a 
value of / is chosen randomly according to the distribution shown in equation 2. Then the 
charges of the pions are interchanged pairwise (tt + tt~ or tt 7i°) until the charge distribution 
matches the chosen value of /. This simulates a DCC accompanied by the normal hadronic 
background in a way that conserves energy, momentum, and charge. The Sz distribution 
for the ( = 0%, 25% and 60% DCC hypotheses are shown in figure 2. The distributions 
get wider as ( is increased. Thus DCC events would appear as outliers with respect to 
the bulk of the data. 

4.2. DCC upper limit 

We expect the DCC events to show up as non-statistical tails in the Sz plot. Since 
we do not see no such events in our data sample, we are faced with the possibilities that 
single-domain DCCs are very rare, very small, or both. To check which hypotheses are 
consistent with our data, we determine upper limits on the frequency of DCC production 
as a function of its size, as represented by (. 

We have calculated the 90% confidence limits for two scenarios and the results are 
shown in the right side of figure 2. The first, shown as the solid line, is based upon the 
conservative assumption that the MC simulation should describe the data perfectly in the 
absence of a DCC signal. The second scenario, shown as the dashed line, assumes that 
the difference between the data and VENUS is due to detector effects and that the widths 



should be the same. Further details can be found in [ [16] 
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Figure 2. (a) Sz distribution for the experimental data is shown, overlaid with simulation 
results incorporating 0%, 25% and 60% DCC in every event, (b) 90% C.L upper limit on 
DCC production per central event as a function of the fraction of DCC pions under two 
assumptions. 



5. Methods for DCC domains localized in phase space 

The localized DCC analyses take advantage of the predictions that the coherent fluc- 
tuations caused by DCC tend to cluster in phase space. A similar analysis technique that 
described in the last section can be used in small phase space. Here we describe two other 
analysis methods, viz., use of wavelets and factorial moments. 

In WA98, the localized analysis is done by taking the data of the full overlap regions of 
PMD and SPMD, which lies between 77 = 2.8 and 3.75. In figure 3 this zone is shown for 
a single event in terms of x — y phase space distribution on the PMD plane. Open circles 
show the distribution of 7-like clusters from PMD and filled circles show the charged 
particle distributions from SPMD projected on to PMD plane. The neutral and charged 
particle distributions are evenly distributed in this event. The methods discussed below 
are sensitive to single out exotic events out of a large data sample. 

5.1. DCC simulation for localized analysis 

For the simulation of normal background events we have adopted the same simulation 
method as described in section 3 whereas the DCC simulation is done differently compared 
the global case. Since the localized methods are sensitive to smaller phase space zones, 
we have the option of introducing single or multiple domains of varying sizes to any phase 
space into the simulation and carry out further analysis. For the present study, we have 
chosen a DCC domain to be within 3 < rj < 4 and < <fi < 90. The common phase space 
of PMD and SPMD overlaps with this choice of domain. The rest of the analysis follows 
as described in section 4.1 including the full GEANT simulation of WA98 detectors. 
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Figure 3. x — y phase space distribution of a normal event from the WA98 experiment on 
Pb+Pb reactions at 158-A GeV/c. 



5.2. Wavelet Analysis 

A unique analysis method based on the discrete wavelet transformation (DWT) is 
adopted to search for the fluctuation in neutral pion fraction in a localized i] — <p space. 
This method was first proposed by Huang et al. [ f° r DCC search. Wavelets are the basis 
functions in some representations of arbitrary functions that satisfy certain requirements 
like invertibility, orthogonality etc. (as sines and cosines in case of Fourier transformation). 
These arbitrary functions can be functional representations of a given data set. The most 
interesting feature of the discrete wavelet transformation is the zooming action at each 
location of the various resolution scale. Due to the completeness and orthogonality of the 
DWT basis, there is no information loss at any scale. 

In our present analysis, we have chosen the arbitrary function to be the neutral pion 
fraction given in equation (1). The DWT representation of the neutral pion fraction can 
be represented by, 

2i-i_l 2^- 1 -l 

f b) (x)= J2 fj-i,k<t>j-i,k( x ) + J2 fj-i,k^j-i,k( x ) ( 3 ) 

k=0 k=0 

where fj^ are the mother function coefficients (MFC) and fj^ are the father function 
coefficients (FFC) at position index k of resolution scale j. The MFC's and FFC's are 
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Figure 4. (a) FFC distribution of simulated data of different percentage of DCC events, 
(b) FFC distribution of data and comparison with MC simulation. 



the carriers of information at each scale. In Haar basis, MFC's are the average between 
two bins and FFC's are the half difference between two bins. The variable x can be r], <fi 
or a combination of both 77 and 0. We have chosen x to be 77 in our case. 4>(r]) and ip{rj) 
are called the mother and father functions, respectively. In the present analysis we have 
used D-4 wavelet basis. 

In order to verify if the FFC distributions are sensitive to events with DCC pions we 
have carried out the MC simulation described above for different percentage of DCC 
events present in the full sample. We show two of these cases in figure 4(a). The solid 
histogram shows the FFC distribution at j = 2 in the absence of any DCC event. This is 
our background simulation. The dotted and dashed curves represent the cases where DCC 
events are present in all the (100%) events and 10% of the events, respectively. In all cases 
the DCC domain is chosen to be the same, as described in previous section. Significant 
difference between the background histogram and events with DCC can be seen here. In 
this method one is sensitive to a small admixture of DCC events in the full sample. It is 
seen that even if the occurrence of DCC events compared to the total number of events is 
very less, still then there are significant number of events beyond 3a of FFC distribution 
of background events at j=l. These events are to be tagged and studied in detail for the 
presence of DCC. 

The WA98 data FFC distribution (at j = 2) for central events is shown as solid points in 
figure 4(b) and overlaid on the top is the background MC simulation. Data and simulation 
match quite well. Detailed analysis is in progress by rotating the full phase space by a 
given angle in <f> and repeating the whole analysis. Another way is to rotate the SPMD 
plane with respect to PMD plane and determining FFC from the resulting distribution. 

Another interesting quantity in the fluctuation analysis at each scale is the wavelet 
power spectrum, defined as: 

Pi = hE\fik\ 2 - (4) 

1 k=0 
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Figure 5. (a) Wavelet power spectrum from simulated MC events for different percentage 
of DCC events (b) Power spectrum for WA98 data. Events beyond the three sigma value 
of background MC simulation are to be selected as DCC candidates for further study. 



Figure 5(a) shows the power spectrum calculation for different cases. The power spectrum 
for background MC simulated events is quite flat at all scales, j. In the left hand panel 
for figure 5 we show simulation for 100% and 10% DCC events. In these cases we are 
averaging over many events which is making the Pj very weak. That means, Pj at j=l 
which is so much higher for 100% DCC events has gone down below noticeable level 
for 10% DCC. Instead of averaging over events which kills the signals this plot will be 
presented event-by-event. 

One can filter out significant number of exotic events by selecting the events which have 
FFC beyond the three sigma level of the fitted gaussian of the solid curve of figure 5(b). 
The power spectrum corresponding to the three sigma value of FFC is shown in the right 
panel of the figure as a solid curve. The selected events will have a power value above 
this line. These events, once found, have to be analyzed further by looking at 77 — or 
x — y event display plots for the existence of DCC. 

5.3. Bivariate Moments 

Since pions from the DCC domains would appear as localized fluctuations in the joint 
charged to neutral particle distributions, one of the elegant analysis methods is to calculate 
bi-variate factorial moments such as: 

f lfi (ch, 7 ) = (N ch ), fo,i(ch,<y) = (N J ), 
f ltl {ch,j) = (N ch N-y), / 2 ,o(c/i,7) = (N ch (N ch - 1)) (5) 

and so on. Here the averaging is done over a large number of samples. Filtering of true 
fluctuations from the above distributions become quite complicated because of various 
detector effects, efficiency and background factors. Some of the deficiencies could be over- 



come by using the generating function technique [ [TSJ] which can be applied to concerned 
distributions. In that way we can evaluate normalized factorial moments for the distribu- 
tions which have the ability to filter out the non-statistical fluctuations from the envelop 
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of statistical noise. One defines the normalized bi-variate factorial moments as: 




(N ch (N ch - 1) . . . (N ch - % + 1) iV 7 (iV 7 - 1) . . . (Ky - j + 1)) 

(N ch ) 1 (N,Y 



(6) 



One defines the ratios of these factorial moments as: 




(?) 



The Minimax collaboration has shown that these ratios are not affected by detector 
effects or efficiencies, and termed as "Robust Observables" . Here the sampling is done 
over total number of events, rather than event-by-event. These ratios yield 1 for generic 
distributions and l/(i+l) for pure DCC. 

In case of WA98 data the these observables are not robust, but one can calculate these 
quantity in case of both simulation with and without DCC and then compare with data. 
In order to make these quantities suitable for identifying single event fluctuations the 
averages are taken over a considerable number of 77 — <p phase space bins of iV 7 and N^. 
This can be studied in WA98 because of the large acceptance of both photon and charged 
particle detectors. In average, the results from the data match with MC simulation. 
Detailed analysis is in progress to sort out events with large fluctuations. 

6. Future Prospects 

First we outline some of the new analysis from WA98 experiment by combining with 
neutral and charged particle spectra. Photon and charged particle momentum spectra 
for different event classes based on iV 7 _iike and N c ^ distributions (e.g., by choosing events 
around the mean value of FFC and events far from the mean) may show significant 
differences since DCC pions are predominantly of low pt- Average of photons obtained 
from PMD may also show differences for these event classes. 

Several other signals of DCC have been discussed in the literature for the detection 
of DCC. It has been proposed that the photon and dilepton signals will be significantly 
enhanced in the presence of DCC. This has been discussed by V. Koch and by J. Randrup 
in this conference. M. Asakawa has argued (also in this conference) that non-central 
collisions may be a better place to look for DCC. The two particle correlation (HBT) 
effect may also be suppressed if a large domain forms [ H. Another interesting signal may 
be the enhancement in the production of baryons and antibaryons [ |19| . 

Most of the future experiments at RHIC and LHC have plans to search for the signals 
of DCC. A highly granular photon multiplicity detector is being planned for the STAR 
experiment at RHIC, which, in combination with charged particle detectors and FTPC, 
will be quite adequate for DCC search. Here one would be able to select on the low 
Pt particles, characteristic of DCC pions. DCC search in PHENIX will be possible by 
correlating signals from charged particle detectors and photons from EMCAL. The PHO- 
BOS experiment, with its large acceptance charged multiplicity detector, will be able to 
search for unusual fluctuations in the production of charged particles. The fluctuations 
can then be correlated to the production of very low pt hadrons measured in its two arm 
spectrometer. The BRAHMS experiment is considering to add a photon arm to look for 
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DCC. A Photon Multiplicity Detector similar to that of STAR is proposed to be included 
in the ALICE experiment at LHC, which along with charged particle detectors will be 
able to search for DCC. 



7. Summary 

The phenomenon of DCC is quite interesting, its observation in heavy ion collisions 
would signal the chiral transition as well as bring a wealth of information for understanding 
of QCD. Since the hope is for such transition to occur at RHIC and LHC energies, we have 
to be well equipped with precision measurement tools and sophisticated analysis methods. 
In this article we have outlined a list of analysis methods that exist at present with a 
speculation on other signals and new techniques for future. With continued development 
and understanding in theoretical aspects and analysis tools, one would hope for the best 
to come from the new data from the upcoming experiments. 



TKN acknowledges very fruitful discussions with Krishna Rajagopal, Jorgen Randrup, 
Ajit Srivastava and Xin-Nian Wang. 
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